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Abstract 

We present experimental evidence that the motion of point defects in thermal con- 
vection patterns in an inclined fluid layer is well-described by Tsallis statistics with 
an entropic index q ~ 1.5. The dynamical properties of the defects (anomalous dif- 
fusion, shape of velocity distributions, power law decay of correlations) are in good 
agreement with typical predictions of nonextensive models, over a range of driving 
parameters. 
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1 Introduction 



In recent years a generalized version of statistical mechanics proposed by C. 
Tsallis, dubbed nonextensive statistical mechanics [1,2,3], has been subject 
to intensive discussion [4]. The basic idea underlying this new theoretical ap- 
proach is quite simple: suppose that a physical system of sufficient complexity 
cannot, for some reason, maximize the usual Boltzmann-Gibbs-Shannon en- 
tropy leading to the usual statistical mechanics. In such a case, it is reasonable 
to conjecture that the system may then maximize some other, more general 
entropy measure. In particular, good candidates for such measures are the 
Tsallis entropies Sq, which depend on a real parameter q, the entropic index. 
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If g 7^ 1, the Tsallis entropies are nonextensive, i.e. nonadditive for indepen- 
dent subsystems. For q = 1, they reduce to the Boltzmann entropy and the 
usual statistical mechanics is recovered. 

If the Tsallis entropies are maximized subject to suitable constraints, one ar- 
rives at power-law generalizations for the canonical ensemble, with q being 
related to the exponent of the power law. An important (and open) question 
involves determining for which types of systems this more general formalism is 
physically relevant. We can imagine various reasons why a particular physical 
system may not be able to maximize the usual Boltzmann entropy: long-ranged 
interactions/correlations, multifractality, metastability, or simply the fact that 
the system is not in equilibrium due to some external forcing. Indeed, recent 
work has shown that nonextensive statistical mechanics is particularly use- 
ful in describing driven nonequilibrium systems. Useful physical applications 
include Eulerian [5,6] and Lagrangian [7,8] fully-developed hydrodynamic tur- 
bulence, heavy ion collisions [9], e+e" annihilation experiments [10,11], as well 
as economic [12] and biological [13] systems. 

Here we report experimental observations of point defect (dislocation) motion 
in defect tTirbiilence in inclined layer convection [14,15], which are found to 
be consistent with Tsallis statistics. To the best of our knowledge, this is 
the first observation of Tsallis statistics in a defect turbulent pattern-forming 
system. The defects behave much like particles obeying a generahzed TsaUis- 
type statistical mechanics, with an entropic index of approximately q ~ 1.5. 
Their dynamical properties (anomalous diffusion, skewness, power law decay 
of correlations) are in good agreement with typical predictions of nonextensive 
models over a range of driving parameters. In addition, since this system is 
but one of many driven systems exhibiting defect turbulence, these findings 
suggest that similar non-extensive behavior may be found in other pattern- 
forming systems such as electroconvection in liquid crystals [16], nonlinear 
optics [17], and auto-catalytic chemical reactions [18]. 



2 Experiment 

We examine the motion of defects within spatiotemporally chaotic patterns 
formed in experiments on inclined layer convection (ILC), a variant of Raleigh- 
Benard convection [19] in which the thin fiuid layer (heated from below and 
cooled from above) is inclined by an angle 9 from the horizontal [14]. For 
intermediate angles of inclination (15° < 9 < 75° at Prandtl number 1), the 
system exhibits defect turbulence within a pattern of undulating convection 
rolls. This state is known as undulation chaos [14,15]; an example is shown 
in Fig. 1, with defects marked by boxes and a tearing region marked by an 
arrow. 
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Fig. 1. Example shadowgraph image of undulation chaos in fluid layer heated from 
below and cooled from above, inclined by an angle 9 = 30°. The nondimensional 
driving paramter is e = 0.08. Black box encloses a positive defect; white, a negative. 
Arrow is adjacent to a tearing region of low-amplitude convection. Uphill direction 
is at left side of page. Region shown is the subregion of size 51d x 63d used for 
analysis. 

The experimental apparatus and data collection for this investigation have 
previously been discussed in [14,15], and only key points are summarized here. 
The planform of the convection pattern was obtained using the shadowgraph 
technique [20] and a digital camera. The fluid used was compressed CO2 with 
Prandtl number Pr = z///t ^ 1, where v is the kinematic viscosity and k is 
the coefficient of thermal expansion. The convection cell had a thickness of 
d — (388 ± 2)/xm and was of dimension lOld x SOd, of which only a central 
bid X 63(i region was used for analysis. The vertical thermal diffusion time 
was = / K = 1.53 sec. All data was collected at an inclination of = 
30°. The strength of the thermal driving is described by the nondimensional 
temperature difference e = — 1; this quantity was varied in the experiments. 
Data collection occurred at 3 frames/sec in two modes: four to six hundred 
100-frame runs at each of 17 values of 0.04 < e < 0.22 and a few (4, 2 
respectively) 80,000-frame runs at e = 0.08 and e = 0.17. Each value of e was 
reached via quasistatic temperature changes. We also conducted a sequence of 
measurements with quasistatic temperature increases followed by quasistatic 
temperature decreases to check for possible hysteresis, which was not observed. 
The Boussinesq parameter [19] was Q — 0.8, indicating a breaking of the up- 
down (z) symmetry. 

Due to the inclination of the fluid layer, the system is anisotropic and the 
convection roll axes align along the direction of the shear flow (y) for interme- 
diate inclinations such as ^ = 30°. In contrast to what is observed in Rayleigh- 
Benard convection, the defect motion is observed to be primarily transverse 
to the convection rolls ("glide," in the language of crystals) instead of parallel 
("climb"). As such, the motion of the defects does not typically adjust the 
local wavenumber of the rolls, rather their orientation and the wavenumber of 
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the undulations. 



3 Nonextensive statistical mechanics of defects 

3.1 Observed defect statistics 

Topological defects in spatiotemporally chaotic patterns such as undulation 
chaos in ILC behave much like particles: they have charge; they are created 
and annihilated; and they have a well-defined position and velocity. Their 
topological charge is a quantity analogous to a Burger's vector; point defects 
(dislocations) are located at positions where there is a discontinuity of ±27r in 
the phase of the pattern $ along a contour around the point: 



As seen in Fig. 1, undulation chaos consists of an underlying set of stripes (con- 
vection rolls) which contain both undulations and such defects. Each of these 
features is observed to be spatiotemporally chaotic. While undulation chaos 
has strong nonequilibrium properties due to the thermal forcing, it nonethe- 
less exhibits stationarity in distributions of the number of defects and the 
wavenumber [14]. 

If the defects behaved like an ideal gas obeying standard statistical mechanics, 
the probability P^r to observe particles in a certain area would be given 
by the Poisson distribution P^r = ^ where A is the average number of 
particles. However, this distribution has been observed to be invalid for a 
number of defect turbulent systems [21,14,22]. A particular model for P/v, 
based on empirically determined rate equations, has been studied in [14]. It 
yields a modified Poisson distribution 



where It, denotes the modified Bessel function and the constants a and b are 
determined by the defect gain and loss rates. This prediction differs signifi- 
cantly in its variance from the Poisson distribution, and agrees very well with 
the experimentally observed distributions in undulation chaos [14]. 

Besides defect number, it is also interesting to look at the statistics of defect 
velocities. We examine the ensemble of defect trajectories decomposed into 
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Fig. 2. Transverse velocity (vx) distributions at (a) e = 0.08 and (b) e = 0.17 for 
positive and negative defects, rescaled to cr = 1. Solid lines are one-parameter fits 
to Eq. (6) for positive defects. Unrescaled standard deviations were: (a) cr+ = 0.550, 
cr„ = 0.553 and (b) cr+ = 0.586, cr_ = 0.565. Inset: Relative error of experiment 

and theory: {Pexp - Ptheory) / Ptheory 
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Fig. 3. Longitudinal velocity (vy) distributions at (a) e = 0.08 and (b) e = 0.17 for 
positive and negative defects, rescaled to cr = 1. Solid line is a fit to Eq. (10) for 
positive defects. Unrescaled standard deviations were: (a) (7+ = 0.099, (T_ = 0.115 
and (b) C7+ = 0.141, c7_ = 0.128. Inset: Relative error of experiment and theory: 

{Pexp Ptheory) / Ptheory 



the two directions x (across rolls) and y (along rolls, uphill-downhill) coin- 
ciding with the anisotropy of the system. This decomposition is done because 
their motions are fundamentally different with respect to the underlying roll 
structure; their dynamics have been observed to evolve in a nearly statistically 
independent way [15]. If the defects behaved like ordinary particles obeying 
ordinary statistical mechanics, one would expect Gaussian velocity distribu- 
tions and normal diffusion of the positions. Instead, significant deviations are 
observed [15]. 
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The probability density functions (PDFs) of the velocities and Vy are shown 
in Fig. 2 and 3, with velocities rescaled by their standard deviation a. Data 
is shown for e = 0.08 and 0.17, the two values with the largest datasets. The 
PDFs deviate significantly from a Gaussian, with velocities of more than 20fT 
observed in the long tails. In the transverse direction (x), the velocities are 
several times faster than those in in the longitudinal (y) direction. 

A common characteristic of nonextensive systems is the presence of long-range 
interactions and/or long-range correlations: for generic nonextensive models, 
velocity correlation functions are expected to decay asymptotically with a 
power law. While the temporal correlation function C{t) of defect velocities 
decays exponentially for small t (< lOrt,) [15], Fig. 4 reveals a superimposed 
long-term power law contribution of the form C{t) ~ t"'' for large t. For both 
Vx and Vy at both values of e, the exponent r) is greater than 1, although the 
fit region is of limited length due to finite observation times. 

Finally, these observed velocities provide for defect diffusion. Normal diffusive 
behavior is given by the second moment of the position {x{tY) — {x{t))'^ oc 
where a — 1 and (•) represents an ensemble-average. Defect motion in 
undulation chaos exhibits super diffusive behavior in both x and y, as shown 
in Fig. 5, where the exponent a is greater than 1. 

3.2 Generalized canonical distributions 

Given this anomalous behavior, we investigate the proposition that the de- 
fects obey a generalized statistical mechanics of the Tsallis type and work out 
the corresponding predictions. Because the number of defects is not constant 
but fiuctuates, such fiuctuations should properly be embedded into a gen- 
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Fig. 4. Tails of velocity autocorrelation functions for Vx {x) and Vy (o) at (a) e = 0.08 
and (b) e = 0.17. Straight lines indicate fits to power laws. These curves are the 
tails of the results shown in Fig. 4 of [15]. 
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e = 0.08 e = 0.17 




Fig. 5. Second moment of position trajectories in x and y. Solid line is diffusive 
behavior predicted by Eq. (18) for q = 1.5; dotted line is normal diffusive limit 
(9 = 1). 

eralized grand canonical description. To date, however, no general theory of 
the nonextensive grand canonical ensemble is known, and we instead use the 
canonical nonextensive model. This theory predicts that the stationary PDF 
of an observable u should be given by the generalized canonical distribution 

p{u)^^{l + {q-l)PE{u)y^K (3) 
Here E{u) is an effective energy associated with u, and 

+00 

Z,^ J {l + {q- l)PE{u))-^^du (4) 

— oo 

is a normalization constant. /5 = l/{kT) is a (formal) inverse temperature 
variable, which essentially fixes the variance of the distributions, and q is the 
entropic index. 

The distribution in Eq. (3) is obtained by maximizing the Tsallis entropies 




subject to suitable constraints [1,2]. The {pi} are the probabilities of the mi- 
crostates i of the system. In the limit q ^ I, the usual Boltzmann distribution 
is recovered from Eq. (3). 



7 



The lowest-order model is E{u) = |it^, which may be associated with an 
effective kinetic energy. For this relation, one obtains the PDF 

p{u)^^(l + hq-l)PuA~~\ (6) 



where 



2 r(i)r(^-|) 



(7) 



The distribution in Eq. (6) has variance 1 if /3 is chosen to be 



yielding a single fit parameter q. 

In general, however, defect interactions, the local pattern, and other higher- 
order correlations of the chaotic forces acting on the defects will produce small 
correction terms to E{u). The effect of such higher-order correlations of chaotic 
driving forces has been studied in [23,24]. The result of those considerations is 
that for distributions which have been rescaled to a = 1 the effective energy 
is of the form 

E{u)^^e + c(u-lu'^+0{c'). (9) 



Here c is a small constant related to the skewness of the distribution and 
whether or not it is nonzero depends on the universality class of the chaotic 
driving forces that are acting. Eq. (9) is a pcrturbative result that is valid for 
u <^ c^^. Neglecting higher-order corrections of O(c^) we obtain the formula 



for the probability distribution. 

When this formalism is applied to defects, we associate u with either or Vy, 
rescaled by {v'^^y}^^^ to obtain an appropriate dimensionless quantity. Concep- 
tually, this corresponds to associating the u'^ term of E{u) with an effective 
"kinetic energy," in the language of ideal gases. In the case of the transverse 
motion (£, Fig. 2) of the defects, the measured PDFs are symmetric and coin- 
cide very well with symmetric (c = 0) Tsallis distributions of the form Eq. (6), 
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Fig. 6. Overlay of (a) Vx PDFs and (b) Vy PDFs for positive defects at 17 values 
of e (open, gray symbols) with data at e = 0.08 (+) and e = 0.17 (x). This data 
corresponds to the peak region of Figs. 2 and 3, but not rescaled by a. 



while for the y motion (Fig. 3) significant asymmetries are present and Eq. (9) 
is appropriate with c 7^ 0. The presence of this skewness in the Vy may be as- 
sociated with the non-Boussinesqness of the fluid flow, which causes the shear 
flow to be asymmetric in z. Previous studies of inclined layer convection [14] 
have revealed similar asymmetries: the undulations drift in the y direction. 

We flt the experimentally obtained PDFs to either Eq. (6) or (10) by min- 
imizing T,{Pexp - PtheoryY /Ptheory, which allows for appropriate weighting of 
the tails. These fits can be seen as the solid lines in Figs. 2 and 3. For v^, 
we obtained values oi q = 1.45 at e = 0.08 and 1.51 at e = 0.17, with the 
value for (3 set via Eq. (8) by rescaling the distributions to o" = 1. For Vy^ we 
find q = 1.40 and c = 0.0125 (such that the expansion is valid for u ^ 80) 
at e = 0.17. For Vy a.t e = 0.08, the data is sufficiently skewed to violate the 
conditions under which the expansion in c is valid, and no good fit was found 
even for values of c > 0.03. It should be noted that these fits do display small 
systematic deviations from the data, which may or may not be significant over 
six orders of magnitude. 

The PDFs for the velocity components are not significantly changing in q, j3, 
and c across the values of e surveyed. This is demonstrated in Fig. 6 for and 
Vy in which the PDFs are found to coincide with no rescahng. Fig. 7 shows the 
fit values for q and c at various e. A mean value of ^ = 1.48 ± 0.05 is observed 
for Vx- For Vy, we observe q = 1.51 ±0.15 and c = 0.012 ±0.002 (with the sign 
of c corresponding to the sign of the defect). Remarkably, the distributions are 
constant in spite of a transition [15] from disordered undulations for e < 0.1 
to intermittently ordered undulations for e > 0.1. 
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Fig. 7. Values of fit parameters (a) q for ■Uj; (b) q for Uj^, and (c) c for i^j; as a function 

of e. Data in black is for positive defects; data in gray is for negative. Short runs 
in which data was obtained by quasistatic temperature increases from below are 
marked with +, and the converse with o. 

3.3 Effective degrees of freedom of the environment 



We now consider concrete dynamical models exhibiting nonextensive behavior. 
One of the simplest models, described in [6], is a linear Langevin equation 
whose friction parameter is also stochastic. We use such an equation as a 
simplified dynamical model for defect velocities: 

ii — —^{x,t)u + KL{t), (11) 



where L{t) is Gaussian white noise of strength k, and 7(0;, t) is a random 
damping constant which fluctuates and may be attributed to spatiotemporal 
changes in the environment of the defect. One can let k, fluctuate as well [6], 
but this is not necessary for our purposes. 

It has been shown in [6] that Eq. (11) generates Tsallis distributions for u if 
P = ^/k^ is x^-distributed with degree n, i.e. if the PDF of /3 is given by 



n(3 
'Wo 



(12) 



Such a distribution is generated by a sum of n Gaussian random variables Xi, 
squared: 

n 

/3 = E^' (13) 



If (3 fluctuates on a time scale that is much larger than the relaxation time to 
reach local equilibrium, then, as shown in [6], Eq. (11) generates stationary 
densities of the form (6) with 
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(14) 



n+V 




(15) 



Here f3o = Jq° f{(3)(3 d(3 is the average value of the fluctuating /3. The ex- 
perimentally observed value g ?^ 1.5 for the defect statistics thus means that 
there are effectively 3 independent degrees of freedom that contribute to the 
fluctuating local defect environment. 

Defects are no ordinary particles: they have neither a well-defined mass nor 
a well-defined size. Furthermore, the nature of defect turbulence means that 
they are moving within a spatiotemporally chaotic environment. Therefore, 
one expects that there is an ensemble of damping constants which depend on 
the environment of the defect, providing the fluctuating effective damping 7. 

In particular, the fastest velocities result from circumstances in which the 
defect is moving in a local environment corresponding to weak damping (small 
7). In undulation chaos, such events coincide with the presence of a transverse 
region of low-amplitude convection, as shown by the black arrow in Fig. 1. The 
rapid motion takes place due to the fact that the pattern "tears" in such a 
way that each roll broken rejoins with the one next to it, rapidly transporting 
the defect across the weak region. The driving forces L{t) are weakly damped 
during such a time interval, leading to very large velocities for short intervals, 
until another region with another 7 is reached. The behavior of such flight-like 
events is further described in [15]. 

In general, one may also consider probability distributions other than the x^- 
distribution in used Eq. (12). These lead to statistics other than Tsallis statis- 
tics: superstatistics [25] , for which generalized versions of statistical mechanics 
can be constructed as well [26]. Using these extensions, it can been shown that 
for /5-fluctuations of small variance and not too large all superstatistics are 
universal, and in fact are given by Tsallis statistics [25]. Significant differences 
arise only for large \u\. In general, the x^-distribution is distinguished by the 
fact that it yields a superstatistics whose generalized entropies have "better" 
properties than those of other superstatistics. 



3.4 Anomalous diffusion 

Typical nonextensive models predict anomalous diffusion, such that 




(16) 
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with a 7^ 1, as is observed in the diffusion of defects in undulation chaos. 
The nonhnear Fokker-Planck model previously discussed in [27,28] describes 
correlated anomalous diffusion in the framework of nonextensive models. It is 
based on a generahzed Fokker-Planck equation of the form 

t) = {F{x)p{x, t)) + D^^{p{x, t)r (17) 



with a linear drift term given by F{x) — ki — k2X. 

The linear case {u = 1) corresponds to the standard Fokker-Planck equation 
generating the usual Brownian motion processes with Gaussian densities. As 
shown in [27], the nonlinear case where i' ^ 1 generates Tsallis distributions, 
with the entropic index q being related to the parameter u hj q = 2 — u. In 
addition, the model generates anomalous diffusion for u ^ 1. with the relation 
between the exponent a of anomalous diffusion and the entropic index q given 
by [27] 

2 2 

= . (18) 



Thus, we can compare the results for the Tsallis fits to the PDFs and the 
exponents in the diffusion graphs. Remarkably, the value a = 1.33 obtained 
from Eq. (18) for q = 1.5 provides a good match to the observed data. Other 
models, based on 2-dimensional Levy flights [29], predict a — q — 1 0.5 and 
are apparently not consistent with our data. This can be understood within 
the framework described here: we earlier found that n ^ 3, which means 
that there are effectively more than 2 degrees of freedom contributing to the 
fluctuating environment of the defect. 



3. 5 Velocity correlation function 



To better understand the relationship between Figs. 4 and 5, we can perform 
a simplified calculation relating the velocities and positions. To do so, we 
neglect all spatial correlations, which do effect the diffusion but have not been 
measured, and focus only on the temporal ones. In each spatial direction the 
displacements obey x{t) = jlv(t')dt'. Hence, 



{x^(t)) = j dt' j dt"{v{t')v{t")) (19) 
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t t' 



2 dt' dT{v{0)v{r)). 



(20) 







Assuming the simple asymptotic form {v{0)v{t)) ~ t the two integrations 
in Eq. (20) yield 



as the expected behavior in the asymptotic hmit. By equating this expres- 
sion with Eq. (16), one obtains a relation between r] (describing the velocity 
autocorrelation) and a (describing the diffusion), namely 



In the experimental data displayed in Fig. 5, we observed that the effective 
anomalous diffusion exponent a is given by a ~ 1.33 on moderate time scales 
and decreases to a < 1 for large time scales t > lOr^. Thus, Eq. (22) yields 
~ 2 — a > 1 for the long-term decay exponent of the velocity correlation 
function. As Fig. 4 shows, this is in good agreement with the measured decay 
rate of the velocity autocorrelation function. 



4 Conclusions 

We have provided experimental evidence that topological defects in Inclined 
Layer Convection behave in a similar fashion to an ideal gas of nonextensive 
statistical mechanics with q 1.5. This is an application to a physical system 
far from equilibrium, where the standard statistical mechanics has little to 
say. Not only do the measured probability distributions of defect velocities 
agree well with Tsallis statistics, but we find quantitative agreement in the 
nonextensive parameter 5 1.5 as determined from both velocity distributions 
and diffusive behavior. 

While the direct observation of nonextensive behaviour of the entropy may 
not be possible , further work may nonetheless be able to provide a prediction 
for q. One promising approach would be to quantify the fluctuations in /3, 
the environment of the defect, in terms of observed pattern properties such 
as local wavedirector and convection amplitude. If so, then the recent super- 
statistics described in [25] could provide a suitable framework. Furthermore, 
other defect-turbulent systems may also exhibit similar behavior and prove a 
fertile ground for further investigations. 





r] — 2 — a. 



(22) 
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